Magnetoroton instabilities and static susceptibilities in higher Landau levels 
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We present analytical results concerning the magnetoroton instability in higher Landau levels 
(LL's) evaluated in the single-mode approximation. The roton gap appears at a finite wave vector, 
which is approximately independent of the LL index n, in agreement with numerical calculations 
in the composite-fermion picture. However, a large maximum in the static susceptibility indicates 
a charge-density modulation with wave- vectors qo(n) ~ l/\/2n + 1, as expected from Hartree-Fock 
predictions. We thus obtain a unified description of the leading charge instabilities in all LL's. 

PACS numbers: 73.21.-b, 73.43.-f, 73.43.Lp 



The complex behavior of two-dimensional electrons in 
a perpendicular magnetic field remains to be fully under- 
stood. If the electrons are restricted to one of the two 
lowest Landau levels (LL's), the system exhibits the frac- 
tional quantum Hall effect (FQHE). This behavior can- 
not be explained in the framework of the Hartree-Fock 
approximation (HFA), which predicts a charge-density- 
wave (CDW) ground statea rather than a homogeneous, 
incompressible liquid of strongly correlated electrons.H 
On the other hand, HFA does appear to capture the 
physical properties in higher LL's as shown by recent 
experimental at filling factors v — 9/2,11/2,..., which 
confirm the predictions that the ground state of this sys- 
tem is a unidirectional CDW ( "stripes" ) .□ 

Here we present a method which unifies the physics 
of electrons in the lowest LL's with that in higher LL's. 
By studying the low-energy collective excitations in the 
single- mode approximation (SMA), we show that the 
magnetoroton minimum, which indicates an instability 
of the incompressible liojjid towards Wigner crystalliza- 
tion in the lowest LL's,ErQ in fact indicates an incipient 
CDW instability in higher LL's. In order to observe this 
feature, it is not sufficient to evaluate the energy disper- 
sion ofpthe collective excitations, as performed in previous 
works,B0 but one must consider the static susceptibility. 
As we will show, the dispersion always presents a mag- 
netoroton minimum at a wave vector q ~ 1 (in units of 
the inverse magnetic length lg), essentially independent 
of the LL index n. However, the static susceptibility 
exhibits a maximum at q ~ 1 for the lowest LL's, a sig- 
nature of the Wigner crystal (WC) instability, whereas 
the maxima in higher LL's arise at qo(n) ~ l/^/2n*4- 1, 
a feature predicted in the HFA for the CDW phase.Q 

The FQHE phenomenon at filling factors v = 1/ (2s + 
1), with-.s being an integer, was first explained by 
Laughlino using trial wave functions. These many-body 
wave functions describe electrons attached to 2s + 1 vor- 
tices, and this attachment leads to a minimization of the 
Coulomb energy, whkjh is the only energy scale in the 
problem. Later, JainS proposed a generalization of the 
trial wave functions, which describe the FQHE at filling 
factors p/(2ps ± 1), where the most prominent plateaus 
are observed experimentally (p is an integer) . The FQHE 
is understood as an integer quantum Hall effect of com- 



posite fermions (CF's) subject to a reduced magnetic field 
and filling p CF-LL's. These new particles are interpreted 
as electrons bound to 2s vortices, and the Laughlin states 
are obtained when the CF's populate the lowest CF-LL. 

If one neglects LL mixing, the underlying Hamiltonian 
describing the low-energy physics of spin-polarized elec- 
trons in any LL n is given in terms only of the Coulomb 
interaction by 



H n =J2v(q) [F n {qf p{-q)p{q), 



(1) 



where v(q) = 2ire 2 /(eq) is the Fourier transfor- 
mation of the Coulomb interaction and F n (q) — 
L n (q 2 /2) exp(— g 2 /4) is the form factor in terms of La- 
guerre polynomials arising from the electron wave func- 
tions in the nth LL. The operators p(q) are the projected 
Fourier components of the electron density, and obey the 
commutation relations 



\p(q),p(k)] = 2? sin 



(q x k) z 



p(q + k). 



LL mixing in the limit where the cyclotron gap HeB/m 
is large compared to the characteristic Coulomb energy 
e 2 /(e?s) may be taken into account in the form of a 
screened interaction potential v(q) = 2ire 2 /e(q)q with 
the dielectric function e(g).a This modification is signif- 
icant only in the vicinity of a characteristic wave vec- 
tor q = 2.4/v / 2rT+T. In the following we will neglect 
this screening because it is unimportanLin determining 
the behavior of electrons in higher LL's.E3 Note that the 
Hamiltonian (|l]) allows one to describe electrons in any 
LL n in terms of electrons in the lowest LL interacting 
via a modified potential v n (q) = v(q)[F n (q)] 2 . 

To describe the collectivepexcitations of electrons in the 
FQHE regime, Girvin et alB used the SMA, in which the 
excited state is given by 



The energy dispersion of the excitations is 
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where E^ s+1 is the energy of the Laughlin state |f2 2s +i)- 
This energy remains finite at any wave vector, as ex- 
pected for an incompressible ground state, but has a 
minimum at a nonzero value of the wave vector. This 
magnetoroton minimum is interpreted as the precursor 
of a WC instability because it matches the reciprocal- 
lattice vector of the WC. The WC is expected to be the 
true ground state at very low filling factors v, where the 
average distance between electrons exceeds the spatial 
extent / ~ lg of the single-particle wave functions. With 
increasing s, the gap at the magnetoroton minimum de- 
creases, indicating that the WC becomes energetically 
more favorable. But in the SMA the gap never vanishes, 
so one would expect the Laughlin state to be the true 
ground state even at very plow filling factors (large s). 
However, Kamilla and Jaintil have shown that the SMA 
overestimates the gap at the magnetoroton minimum and 
fails at larger wave vectors. They compare the SMA to 
a description of the collective excitations as-Jiiagnetoex- 
citons of CF's. The excited state is given bj£3 

where |f2 p g ) is the CF ground state, in which CF's con- 
taining 2s vortices populate p CF-LL's, and pZ~^ p+1 cre- 
ates a particle-hole pair of wave vector q, with the particle 
situated in the (j> + l)th CF-LL and the hole in the pth. 
V is the projection operator to the lowest LL. Even if it 
is not evident how to relate the phononlike excitations 
introduced in the SMA to the CF magneto-excitons_Jthe 
two have been shown to coincide in the limit q —> 0.t3 In 
the CF picture one finds that for v = 1/9 and lower fill- 
ing factors the collective excitations at the magnetoroton 
minimum have a lower energy than the Laughlin state. 
This indicates that the WC is-jthe true ground state of 
the system in the dilute limit Jid 

The question of whether the FQHE states in higher 
LL's become unstable towards a CDW was discussed by 
Scarola et ai.0 using similar arguments as for the WC 
instability. The authors found that the CF's with two 
vortices ( 2 CF's) are unstable for n > 2, because the 
CF excitons have a lower energy. However, the insta- 
bility is still found at wave vectors of order 1, and not 
at the characteristic wave- vectors qo (n) ~ 1/ \J2n + 1 as 
expected for the CDW. Their second remarkable result 
was that excitons of CF's with 4 or 6 vortices ( 4 CF's or 
6 CF's) show a finite gap at all wave vectors for the LL's 
n = 2 and n = 3. This supports the conjecture that 
the FQHE could be found in ultrahigh mobility samples 
around v = 1/4 or v = 1/6 in higher LLIs, where v now 
denotes the filling factor of the last LL.&3 

Here, we present analytical results for the SMA in 
higher LL's. Studies for n = 1 were fixst performed nu- 
merically by MacDonald and Girvin.tl The phononlike 
collective excitations in a higher LL are described by the 
wave function 

\i>nt) = F n {q)p{q)\n 2a+l ), 



where we have made use of the projection scheme pre- 
sented above. One finds for the energy dispersion 

A&Uff) = 2^«(fc)[F„(fc)] 2 sin^^^^ 

x e l3- +g l 2 / 2 s 2s+1 (g + | - e fc2 / 2 s 2s+1 (fc) 
e9 2 /2 S2s+l(g ) 



where S2 S +i(s) is the projected static structure factor of 
the (2s + 1)-Laughlin state. The structure factor is the 
Fourier transformation of the pair distribution function, 
which may be calculated using Monte Carle, integration^ 
The latter may be expressed as the seriesLLJ 

S2.+i(r) = l-e- r2 l 2 

oo _ / 2 \ 2m+l 

+ V c 2s+1 I — I e- r2/A 

+ ^(2m+l)! 2m+1 U/ 

which leads to the formula 

m=0 

for the projected static structure factor. In contrast to 
Ref. 5, where the expansion coefficients c 2 ^^ were ob- 
tained by a fit to the Monte Carlo results, here we use 
the following sum rules as determining equations for the 
coefficients, 

oo 

E 2s+1 _ _£ 
L 2m+1 — 2' 

m=0 

oo 

^(2m + 2X+ + 1 1 = -|, 

m=0 

oo 2 

^(2m + 3)(2m + 2) C ^+ + 1 1 = i-, 

c 2m+l = ~ 1 f° r m < s - 

Coefficients with an index m > M = 3 + s, which is 
the maximal index one can determine in this procedure, 
are set to zero. The first three sum rules are due to 
charge neutrality, perfect screening, and a compressibility 
sum rule, respectively. The last condition is given by 
the repulsion of— the electrons at short distances, where 
g(r) ~ r 2 ( 2s+1 ).EJ The results obtained in the lowest LL 
for pair distribution function, projected structure factor, 
and energy dispersion (see curves in crosses in Fig. |l|) are 
in good agreement with Monte Carlo calculations. 

In higher LL's (see Fig. Q), we obtain the following re- 
sults: for v = 1/3 (s = 1), the gap at the roton minimum 
is largest in the lowest LL (n = 0), n = 1 and n = 2 
determine an intermediate region (not shown) where the 
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FIG. 1: Energy dispersions for v = 1/3 (a) and v = 1/5 (b). 



SMA breaks down, and for n > 3, the roton gap de- 
creases continuously. At v = 1/5 (s — 2), the roton gap 
increases from n = to n = 1, the intermediate region 
is found at n = 3 and n = 4, and for n > 5, we observe 
the decrease of the gap. The same qualitative behavior 
is found for v = 1/7 (s = 3) with a maximal roton gap 
in the n = 3 LL, the breakdown of the SMA at n = 4 
and n = 5, and the final decrease for n > 6 (figure not 
shown). For P =1/5 and P = 1/7, the roton gap first in- 
creases with increasing n [see the n = and n = 1 curves 
in Fig. 0(b)] due to the fact that the system is further 
away from the WC regime. The WC transition appears 
at Dwc{n) ~ l/(2n + 1) because the average distance 
between electrons must be considerably larger than the 
size of the wave functions, which scales as \J2n + 1 in 
higher LL's. Therefore, the WC instability indicated by 
the roton gap is less pronounced. 

For v = 1/5, we find a second minimum in the disper- 
sion relations at q ~ 2.5, see Fig. r-Ub). This minimum 
is already present in the lowest LL.LL3 However, the SMA 
is known to be less reliable at larger wave vectors, and 
for this reason no results were shown in Ref. 5 at q > 2.0 
for v — 1/5. The second minimum also arises in numer- 
ical calculations within the CF picture, which is more 
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FIG. 2: Static susceptibilities for v = 1/3 (a), v — 1/5 (b), 
and v = 1/7 (c). 



accurate in the large-g limit 

Note that we do not find the signature of a CDW state 
in the dispersion relation of the collective excitations, be- 
cause the magnetoroton minimum occurs at wave- vectors 
q ~ 1, independent of the LL. This is-jn agreement with 
the results found in the CF picture.Q The fact that we 
obtain the most stablei-p =1/7 FQHE state in the n = 3 
LL, and not for n — 2,0 is probably due to the proximity 
to the intermediate region in which the SMA becomes 
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less reliable. Breakdown of the SMA in the intermediate 
region is expected to be due to a numerical effect, the 
origin of which is as follows: the gap at q = 



(2) 



may be calculated analytically using the integrals 

II = / dq{k(k-l)L k _ 2 (q 2 )+k(q 2 -2k)L k _ 1 (q 2 ) 
Jo 



^- - (k + l)q 2 + k(k + I) 



L k {q 2 )\[F n {q)] 2 e-^ 



Inspection of /Jm+i as a function of m reveals that the 
integrals are initially negative, and decrease until a min- 
imum is reached around n ~ m. iJm+i then jumps 
to large positive values before decreasing again towards 
zero. It is therefore in the region n ~ m that the gap 
at q = [Eq. (g)] is most sensitive to approximations 
made in determining the coefficients c 2 ^i_ 1 , either nu- 
merically or in the present approach. Because we have 
set Cj^Il! = for m > M, it is in the intermediate region 
defined by n ~ M that the results of the SMA become 
unreliable even in the small-g limit. 

However, the essential aspect of this work is that even if 
one does not find any signature of the CDW in analyzing 
the position of the magnetoroton minimum, it can in fact 
be seen in the static susceptibility calculated in the SMA, 



„ 9 fL+l(g) _ s 2s+1 (q)[L n (q 2 /2)] 2 



A"(g) 



where S2 S+ i(q) is the static structure factor projected to 
the LL n, Girvin et ai.cl have already emphasized that 



the maximum arising in the static susceptibility also in- 
dicates the WC instability in the lowest LL. Because in 
higher LL's the expression for X^s+iil) * s essentially the 
same as for the lowest LL multiplied by the Laguerre 
polynomials of the form factor, the position of the maxi- 
mum now scales as qo(n) ~ l/y/2n + 1 as expected for a 
CDW instability (see arrows in Fig. |^) . 

Note that this interpretation of the CDW in higher 
LL's as arising from the overlap between wave functions 
described by the form factor F n (q) is in agreement with 
the ground-state calculations in the HFA.Q At the wave 
vectors where the Laguerre polynomials have nodes, the 
repulsive Hartree potential vanishes and only the attrac- 
tive exchange interaction remains finite. This gives rise 
to an effective attraction between electrons in the top- 
most LL at the wave- vector go (n) of the CDW. 

In conclusion, we have found a tendency of the FQHE 
states described by Laughlin's wave functions to form a 
CDW in higher LL's. The analytical calculations were 
performed within the SMA. The tendency toward CDW 
formation cannot be observed if one investigates the mag- 
netoroton gap in the dispersion of phononlike collective 
excitations, but becomes apparent only in the static sus- 
ceptibility, which exhibits maxima at the expected CDW 
wave- vectors qo(n) ~ l/y/2n + 1. Our results thus con- 
firm the prediction of CDW formation in higher LL's in 
the framework oiu HFA and provide a unified description 
for the WC and CDW instabilities of all LL's. 
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